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Abstract 

fvg We present a microscopic computation of the radio emission from air showers initiated 

^-H by uhra-high energy cosmic rays in the atmosphere. The strategy adopted is to compute 

P^ each secondary particle contribution of the electromagnetic component and to construct 



the total signal at any location. SELFAS2 is a code which doesn't rely on air shower gen- 
erators like AIRES or CORSIKA and it is based on the concept of air shower universality 
^ which makes it completely autonomous. Each positron and electron of the air shower is 

cn generated randomly following relevant distributions and tracking them along their travel 

^ in the atmosphere. We confirm in this paper earlier results that the radio emission is 

mainly due to the time derivative of the transverse current and the time derivative of 
the charge excess. The time derivative of the transverse current created by systematic 
deviations of charges in the geomagnetic field is usually dominant compared to the charge 



'"5 excess contribution except for the case of an air shower parallel to the geomagnetic field. 
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1. Introduction 



In 1965, Jelley et al. [Ij and Allan in 1971 [2j established that extensive air showers 
00 (EAS) induced by ultra high energy cosmic ray emit a brief radio emission in the range of 

M" few hundreds of MHz during its development in the atmosphere. In the last years, progress 

1/-J made in electronics have permitted the realization of extremely fast radio-detectors used 

^ in the CODALEMA [3], LOPES |1| and AERA |5] experiments. It is clear since ^UU\ 

O and recently P, [3] that the Lorentz force induced by the geomagnetic field acting on 

positrons and electrons of the shower is the dominant source of the EAS radio emission. 
The consequences of this effect are visible through a strong asymmetry in the counting 
rate as a function of arrival directions (see for instance [3l IH |9]). Recent experimental 
analysis shows that the lateral profile of the electric field deposited by an EAS on the 
C^ ground permits to give a good estimation of the primary energy [TU| ITT] . 

In parallel to these experimental results, there are also various theoretical approaches, 
from simple fast semi-analytical models (see Scholten et al. [I2], Chauvin et al. jl3]) to 
more detailed models (see Huege et al [S]). Initially based on the synchrotron radiation 
proposed by Falcke and Gorham [15], the electric field in REAS is now computed using 
the end point formalism [16] applied on each secondary electrons and positrons of the 
air shower. The total electric field at a given observation point is obtained summing all 
electric field contributions coming from the shower. The geometry of the shower (pancake 
thickness, lateral extension) is then taken into account. This characteristic appears to be 
important if we want to compute the electric field close to the shower axis. In MGMR 
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[T7t [12] (Scholten et ai), the total electric field is sourced by the global macroscopic 
transverse current due to the systematic deviation of electrons and positrons in the geo- 
magnetic field. Maxwell equations are applied to this global current concentrated along 
the air shower axis, taking into consideration the thickness of the air shower pancake. 
Due to the secondary electrons excess during the air shower development (mainly due to 
positrons annihilation and knock-out electrons during the air shower development in the 
atmosphere), the residual negative charge variation induces a second contribution to the 
total EAS radio emission [T8| [T9| [IH [20] . This second contribution (charge excess contri- 
bution) is not EAS arrival direction dependent while it is the case for the contribution due 
to systematic deviation of particles in the geomagnetic field (geomagnetic contribution). 
This characteristics implies that the fraction of charge excess contribution to the total 
EAS radio signal also depends on the arrival direction. In most cases, the geomagnetic 
contribution is dominant, but for arrival directions close to the magnetic field orientation, 
the Lorentz force due to the geomagnetic field vanishes and the charge excess contribution 
can become dominant in the total EAS radio signal. 

Our approach is based on a microscopic description of the shower using the concept 
of "age" and "shower universality" first proposed in [21] to study the longitudinal devel- 
opment of purely electromagnetic showers. The use of the relevant distributions for EAS 
secondary electrons and positrons extracted from [221 l23l |2^ (longitudinal profile, particle 
energy, vertical and horizontal momentum angle, lateral distance, and time distribution 
of the shower front), permits to avoid the heavy use of EAS generators to generate air 
showers in SELFAS2 and makes the simulation completely autonomous. Thanks to these 
distributions implemented in SELFAS2, no large simplifications are made on the charac- 
teristics of the electromagnetic shower component. The particles generated in SELFAS2 
by Monte Carlo simulation, are tracked along their trajectory to compute their individual 
electric field contribution to the total electric field emitted by the air shower. With this 
approach, the characteristics of the evolutive spatial density of charge (emissive area) in 
the shower and the systematic drift of electrons and positrons due to the geomagnetic 
field are then naturally taken into account in SELFAS2. 

In section 2, we develop the formalism adopted in SELFAS2 to compute the electric 
field emitted by a single charge with a finite life time which undergoes an acceleration 
due to the presence of a magnetic field. We will describe the strategy adopted to sum 
correctly the contributions given by each particle during its existence in the air shower. 
In section 3 we describe how the electromagnetic component EAS is generated to give to 
each particle initial conditions for their travel in the atmosphere. Section 4 is devoted to 
a discussion from a first example of a vertical 10^^ eV air shower. 

2. Theory: Maxwell equations for moving charges 

2.1. Flash back on modelings 

Historically, many modelings of EAS radio emission were based on the well known 
equation of the electric field emitted by a relativistic moving charge which undergoes an 
acceleration. This equation is obtained solving the Maxwell equation for a retarded time 
using retarded potential (see for example [25j). The resulting electric field of a charge 
moving with a velocity f3 in the lab frame is given by: 
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which is obtained at the observation point x at the time t = tret + R{tret)/c with n = 
{x — r{tret))/R{tret) cind R = \x — r(tret)\- Two temis appear in this equation which 
correspond to a coulombian contribution for the first one and to a radiative contribution 
for the second one due to the acceleration /3 of the source. 

In many approaches ( |26l [27t [28t [T3| |29J), the total electric field generated by the 
EAS is obtained after summation of the individual electric field performed using EqjT] or 
simply replacing q in Eqjllby the global varying charge Q{t) of the EAS. But, as it was 
recently discussed in |20j this way to describe the present problem is not correct because 
EqjT] is obtained for a charge q which is not time dependent. To overcome this problem, a 
solution has been proposed in |14t [T6] . The solution proposed in this paper and adopted 
in SELFAS2 is different. 



2.2. Electric field of a point like source with a finite life time 

Starting from Maxwell equations and using the wave equation for the scalar potential 
$ and the vector potential A (see for instance [25j) we have: 

V^'l'-— ^iTT^ = -— (2) 
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where p is the charge and J the current densities. Using the expression of the electric 
field: 
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we can express the electric field £^ as a function of the charge p and the current densities 



J: 
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This differential equation can be solved using a retarded solution (Green function) 
gives: 
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where we can define R = x — x' and R = \R\ with x' the source position at the retarded 
time t' and x the obervation point. Here, V' and ^ must be considered at retarded time. 
This remark is important, because it is also the case for the expression of the field in Eqj4] 
which should be considered at the retarded time; the variable t in Eqj4] is the time at the 
instant of emission. A particular treatment must be done in order to get V and ^ out 
from the retarded brackets in order to obtain a final expression of ^ as a function of t 
and not t'. Due to the dependance of t' on x' given by the relation between t' and t: 
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the expression [Vp]j,g^ can be transformed in (see [25j 
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where n = R/R is the unit vector between the observation position and the source, 
oriented toward the observation position. Then Eq|6] becomes: 
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Using the fact that the charge distribution is spatially localized on a point, the first term 
in this expression can be rewritten: 



j d'x'^V [p{x', t')U = - j d?x'^^ [p{x\ t')] 
We can finally express the electric field as: 
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which is known as the Jefimenko's generalization of the coulomb law (see [25] fore more 
details). The idea is now to express this formula in the case of a moving particle with a 
finite life time. The expressions for the charge and the current are: 
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where the retarded instant ti corresponds to the creation of the moving charge (particle) 
by sudden acceleration from the state of rest to v and where the retarded instant ^2 
corresponds to the cancellation of the charge by sudden deceleration from v to the state 
of rest. 

Injecting EqjT2]and EqjTSjin EqjTTJand using the fact that R doesn't explicitly depends 
on t we perform the integrations over space and time. We obtain: 
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q{tret) = q[d{tret " ^l) - 9{tret " ^2)] 
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The electric field emitted by a moving charge with finite life time is the summation 
of three contributions. The first one is a static contribution directly linked to the scalar 
potential. It appears by the simple fact that the particle exists. The second term, which is 
probably the less evident term to guess when we see Eqj4| comes from the transformation 
of the retarded scalar potential gradient. It corresponds to the direct charge variation 
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contribution. The third one is directly the time derivative of the current. For the case of 
an isolated particle, it corresponds to its radiation emitted if its undergoes an acceleration, 
like a deflection for instance. For a particle with a finite lifetime, as described by EqfT5| 
the charge variation appears only at the start point and the end point of the particle 
trajectory. Due to the description of the charge existence with a Heaviside-step function 
in Eq|T5| the instantaneous charge variation at ti and ^2 will give un-physical divergences 
for the calculation of the instantaneous electric field at these instants. However, these 
un-physical sharp pulses at ti and t2 can be ignored due to the fact that they give electric 
field contributions out of the frequency range of interest (in our problem, we focus on 
electric field below 500 MHz). We will see in the next section that for an ensemble of 
charge, the collective effect will give a charge variation which, now, has a contribution in 
our frequency range. 

2.3. Electric field of an ensemble of charges 

In the air shower, particles are permanently created and annihilated following the 
distribution of particles number during the air shower development. We compute the 
contribution of each particle during its life time in the shower and we do the summation 
at a given observation point. 

The total electric field created by an ensemble of charges observed at the position x 
is given by: 
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where ( is the number of particle with r{tret) and tret verifying t = tret + \x — r{tret)\/c. 
Performing the time derivative of each particle contribution in EqJT4] before doing the 
summation is not the good strategy. The formulation of Eq|l4] permits us, by a simple 
law of derivation composition, to move the summation symbol inside the time derivative 
operator: 
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This equation summarizes the strategy adopted in SELFAS to compute the electric 
field at any observation point. The strategy will be to construct separately the three 
parts of this equation and to realize the time derivative operation after all particles have 
been considered. The time derivative of charge makes now perfect sense in our frequency 
domain due to the collective effect. The macroscopic charge variation induced by the 
evolution of the number of particles along the air shower development will give a low 
frequency contribution to the electric field. 

To summarize, we can finally decompose the EAS radio emission process into various 
contributions: 

• the static contribution of all particles (term St in the following equation); 

• the time derivative of the net charge excess (term C in the following equation); 
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the time derivative of the transverse current due to systematic separation in the 
geomagnetic field of electrons and positrons from the charge excess and from the 
charge symmetry: electrons and positrons in equal quantity in opposition to charge 
excess (term Cu in the following equation). 



For convenience, we define: 



sti{x,t) 



riiqiit 



ret I 



i?f(l-/3,.n,)J 



and St{x,t) = y ^ sti{x,t) 



(18) 



ret 



{X,t) 



riiqiit 



ret) 



[Ri{l-(3i.ni)\ 



and C{x,t) = y^Ci{x,t) 



ret 



(19) 



CUi{x,t) 



ViQiit 



ret) 



_Ri{l - iSi.ni 



and Cu{x,t) = ^ ^cUi{x,t) 



(20) 



j=i 



Where st{x,t), c{x,t), cu{x,t) are the contributions of the i^^ particle to St{x,t), 
C{x,t), Cu{x,t) which are the total static contribution, the total charge contribution 
and the total current contribution. The total field is then: 



Etot{x,t) = -^ (st{x,t) + ^^^C{x,t) - ^^^Cu{x,t) 



(21) 



3. SELFAS2 algorithm, air shower generation 

3.1. Air shower generation, longitudinal profile, universality 

What we need is to track each particle generated in SELFAS2 along its trajectory to 
compute the electric field as a function of time. To save a lot of computation time, only a 
fraction of the real particle number is generated. The statistical weight and the trajectory 
length are the same in SELFAS2 for each particle generated. The statistical weight is 
calculated by a simple proportionality law between the particle number generated in the 
simulation, and the real particle number which should be generated to obtain the complete 
air shower requested. To respect the distribution of the instantaneous number of particles 
in the shower N{X) (where X is the depth crossed by the shower), we can easily show 
that the random injection depth given to each particle can be performed by Monte Carlo 
simulation following : 
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I 



(22) 



where /, the length of the individual particle track, is sufficiently small with respect to 
the shower length. 

In SELFAS2, the evolution of the particle number is given by the Greisen-Iljina-Linsley 
(GIL) parameterization [3U] based on the Greisen relations [31] and on a variant for de- 
scribing nucleus-initiated showers [7j . This parameterization gives the number of electrons 
and positrons as a function of the primary energy (Ep) and its nature (or mass A): 
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Figure 1: Left: Energy distributions of electrons and positrons in EAS for different evolution stages. 
Right: Charge excess characterization as a function of particle energies for different evolution stages of 
the EAS. The complete set of paranieterizations for particle energy distributions, vertical and horizontal 
momentum angle distribution, lateral distance distributions, and time distributions of the shower front 
is available in 1241. 
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with El = 1450 MeV is a norniahzation factor, E^ = 81 MeV is the critical energy from 
Greisen formula, a = 1.7 is an offset parameter, b = 0.76 is the value of the elongation 
rate coming from adjusted data. X here is the atmospheric depth measured from the first 
interaction point Xi in gcm~^ and Xq = 36.7 gcm~^ is the radiation length of electrons in 
air. The first interaction length Xi can be fixed in SELFAS2 or can be chosen randomly 
following the cross section of the primary with air [32j. 

Using these parameterizations, the evolution of the number of electrons and positrons 
is then completely integrated in SELFAS2, and is completely defined with the energy and 
the nature of the primary fixed in the input file by the user. 

3.2. Particles initial conditions 

With a large set of EAS simulated with CORSIKA, Lafebre et al [24j provide multi- 
dimensional parameterizations for the electron-positron distributions in terms of particle 
energy, vertical and horizontal momentum angle, lateral distance, and time distribution 
of the shower front. They parameterize all these distributions as a function of the relative 
evolution time, t (defined differently as in previous section for the GIL parameterization): 



t 



X — X„ 



(25) 



where X is the integrated thickness of the atmosphere crossed, Xq = 36.7 gcm~^ is 
the radiation length of electrons in air and Xmax is the atmospheric depth for which the 
number of particles in the air shower is the maximum. As an example, the evolution 
of the energy distribution for electrons and positrons is shown in FigjT] left for different 
relative evolution stages. A parameterization of the charge excess is also available in 
|24] as a function of the particle energy and the relative evolution stages of the EAS. 
We show some results in FigjT] right, for different evolution stages. We see that the 

7 



charge excess decreases significantly at high energy (above 100 MeV) during the EAS 
development. When an electron is generated in SELFAS2 with a random energy e, at a 
random evolution stage t, the corresponding charge excess fraction is calculated. Then, 
with a simple random law, the electron is considered coming from the charge excess or 
not. 

All the geometrical parameterizations described in [24j are implemented in SELFAS2. 
At each particle generation, SELFAS2 gives randomly a set of initial conditions following 
all these distributions. The universality of the analytical geometrical descriptions permits 
to generate air showers in the range of 10^^ to 10^° eV with a good accuracy for showers 
coming with zenith angles from 0° to 60°. Only the general characteristics of the required 
air shower is needed by SELFAS2: primary energy, azimuthal and zenith angles. The 
geographical characteristics of the site where the simulation must be realized have to be 
specified as well: orientation of the geomagnetic field, its intensity and the ground altitude 
(the altitude zero corresponds to the sea level). 

3. 3. Particles propagation and electric field computation 

In order to compute energy losses and multiple scattering of electrons and positrons 
in the atmosphere, the length /, of an individual particle track is divided in short tracks 
corresponding to a depth fixed to X^^°^^ (the choice of the resolution is discussed in the 
next section 4). Between two consecutive points, each particle undergoes a deviation 
due to the geomagnetic field and multiple scattering. These two processes are computed 
independently and the final new position in space after one short track results from the 
combination of the two deviations. The geomagnetic field deviation is calculated solving 
the equation of motion and the deviation from multiple scattering is performed by Monte- 
Carlo simulation following the theory of multiple scattering through small angles |33] . 

To compute the electric field contribution along each particle trajectory, the position 
and the velocity of the particle are calculated every X^^°^-^, so we are able to compute 
the instantaneous sti{x,t), Ci{x,t), and cui{x,t) along the particle trajectory, received 
at the time t = tret + R/c at the observation point. On FigM the retarded times after 
each X^^°^^ step for a given particle are represented by tl^^, tl^^ ,tf^^, t^^^ ... To construct 
the shower signal at a given observation point x (an antenna position fixed by the user), 
each particle contribution is added into three histograms corresponding to the three total 
contributions St, C and Cu as a function of reception time t. For these three histograms, 
the constant time width A^'^'l of bins is fixed by the user (tipically 0.1 ns < A^™ < 1 ns). 
For a given particle, to each instant t^^^ is associated a reception time P = tl,^f. + W {tl,^^) / c 
; to fill bins of the three histograms, different configurations are possible (we present the 
method only for the case of Cu histogram but this method is the same to fill St and C 
histograms) : 

• for bins fully-contained between P = t^^^ + R^ (tl^^) / c and P^^ = t^^j +i?-'+-'^(t^^j )/c, 
bins are incremented by : 

CUi{x, P) + CUi{x, P+^) 

e * [Zo) 



with e = 1. Eqj26] means that in SELFAS2, the electric field emitted by a particle is 
considered constant over time during a short track {X^^°^^). This approximation can 
be made because we consider here that the distance R to the observation point, is 
large compared to the size of the particle trajectory during a short track (Fraunhofer 
approximation) . 
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Figure 2: The length of an individual particle track of the particle trajectory is divided in short 
tracks allowing us to compute energy losses and deviations due to geomagnetic field and to multiple 
scattering. At the starting point and at each end of short track, sti{x,t), Ci{x,t) and cui{x,t) are 
calculated independently and added to three independent histograms corresponding to St{x,t), C{x,t) 
and Cu{x,t). Once all particles have been considered, the time derivatives for C{x,t) and Cu{x,t) 
are performed and the complete electric field at this position is obtained summing up the three total 
contributions. 
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• if P is contained into a bin b, delimited by the time t{^;^ and t^^^ + A^™ and if 
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and finally, if P^^ is contained into a bin 6, delimited by the time t{^jj^ and tj^j^^ + A^', 



and if P < t^^^, then the bin b is incremented by Eq,26 with a factor e equal to : 
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Once all particle trajectories have been considered, the numerical time derivatives for 
C{x,t) and Cu{x,t) are performed and the complete electric field corresponding to this 
observation point is obtained summing up the three total contributions. 



4. Discussion from a first example: a vertical 10^^ air shower 

To compare the predicted observables from SELFAS2, REAS3 and MGMR (at the 
end of this section), we simulate a vertical 10^^ eV induced air shower in the Auger site 
configuration. The geomagnetic field characteristics at the Auger site are \B\ =23 /xT, 9b 
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Figure 3: Left: Signal in the EW polarization for an observer located at a distance of 100 m north of 
the shower core. Right: amplitude of the signal as a function of the length of individual particle track 
and its resolution (short track length) fixed in SELFAS2. 



= 58° and (pB = 0° where 6b and (ps are the zenith angle and the azimuthal angle of the 
geomagnetic field. The geomagnetic field comes from the south and is oriented upward. 
The ground shower core of the air shower simulated is located exactly at the center of a 
dense array composed of 145 antennas. Antenna positions are shown in Fig|8]left. The 
ground altitude is 1400 m. 

4.I. Numerical stability 

The final result of the observed signal must be independent of the choice of the particle 
trajectory length discussed in the previous section. With the configuration described 
above, we show in Figj3] left the signal for an observer located at a distance of 100 m 
north of shower core. In Fig|3] we show the amplitude of this signal as a function of the 
length I of the individual particle track and its resolution (short track length) fixed in 
SELFAS2. We observe that for small individual track length (smaller than 5 gcm~^), the 
amplitude of the signal is not stable. In |24j the angular distributions of particles are 
given for both electrons and positrons without discern their nature ; to model the effect 
of the geomagnetic field, electrons and positrons must be propagated. Figj3] shows that a 
minimum length for the individual particle track (/ >10 gcm~^) is necessary to solve the 
absence of geomagnetic field effect in |23] and to guarantee the stability of the signal. 

The choice of the length of the individual particle track is then made in order to obtain 
stability but also to stay as close as possible to the parameterized distributions given in 
|24] . In SELFAS2 this length is fixed to 15 g.cm~^ with a resolution of 0.3 g.cm~^. 

4-2. Results 

In Figjl] we show the east-west polarization radio pulses in the time domain for a 
vertical shower with a primary energy of 10^^ eV computed with SELFAS2. The ground 
frame is arranged so that the air shower core is located at x = and y = 0. The four 
figures represent the pulses observed by antennas set along axis oriented toward the west, 
the east, the north and the south. In these figures we clearly see that the decay of the 
pulse amplitude along the axis west-east is not symmetric around the air shower axis. This 
effect is more visible in Figj5] left where we show the lateral distribution of the absolute 
field strength pulse amplitude. This asymmetry directly comes from the infiuence of the 
net charge excess (more details in the next subsection). The spectral composition of the 
pulses observed in Fig|4] is obtained by Fourier analysis. The result is shown in Figj5] 
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Figure 4: East-west polarization radio pulses for a vertical shower with a primary energy of 10^^ eV 
computed with SELFAS2. In each figure, we show the results obtained for different distances to shower 
core (100 m, 200 m, 400 m and 800 m). Top, left: antennas located along the semi-axis starting from the 
ground shower core pointing toward the west direction. Top, right: antennas located along the semi-axis 
starting from the ground shower core pointing toward the east direction. Bottom, left: antennas located 
along the semi-axis starting from the ground shower core pointing toward the north direction. Bottom, 
right: antennas located along the semi-axis starting from the ground shower core toward pointing the 
south direction. 
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Figure 5: Left: lateral distribution of the absolute field strength pulse amplitude for a 10^^ eV vertical 
air shower. Right: frequency spectra of a 10^^ eV vertical air shower seen from different distances to the 
shower core. East-west polarization. 



right. We see in this figure that the frequency where the coherence is lost, decreases with 
distance to the shower axis. 

To illustrate the notion of coherence loss, let's compute the evolution of energy de- 
posited by the air shower as a function of the number of particles injected in the simulation. 
Every five million particles injected, we performed the spectra integrations in the range 
1-10 MHz and in the range 700-900 MHz (see FigJ6]left). In Fig J6] right, we show the 
energy deposited behavior as a function of the number of particles for the two frequency 
ranges and for different distances to air shower axis. We see that the energy deposited in 
the low frequency range depends on the square of particles number whereas the energy 
deposited in the range 700-900 MHz depends on the number of particles (see comparison 
with 1/A^ in FigJ6] right). This phenomena, due to the coherence loss, is naturally taken 
into consideration in SELFAS2. 



vertical, lO^'eV, SELFAS2 



Coherence vs coherence lost 




500 1000 



1.2 
1.0 

S^O.S 



1/N 



^-■Jj#»..«.j r j-J-4.J-J-J>A.^~*^~I 




Distance / Fcv range 
lOOmll-lOMHz) 
200mh-10MHz) 
400ma-10MHz; 
SOOmh-lOMHzI 
100 m i700-900MHz) 
200 m (700-900MHZ) 
400 m (700-900MHz( 
800 m (700-900MHZ) 



20 40 60 

N Particle number injected [Million] 



Figure 6: Left: Spectra used for the calculation of the energy deposited as a function of the number of 
particles injected in the simulation. The integrations of spectra are performed every five million particles 
injected, for the two frequency ranges 1-10 MHz and 700-900 MHz. Right: energy deposited as a function 
of number of particles, for different distances to air shower axis. The difference of behavior between the 
two frequency ranges of integration is due to the coherence loss of the signal from a particular frequency 
which decreases with distance to air shower axis. 
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4-3. Transverse current vs charge excess 

In section 2, we showed that the EAS radio emission is the summation of three terms. 
The static contribution, St{x,t), is approximatively two orders of magnitude lower than 
the two other contributions, so in the following, this contribution will be not discussed. 
Finally, the time derivative of the transverse current and the time derivative of the net 
charge excess are in competition. In most cases, the transverse current contribution 
is dominant with respect to the net charge excess contribution. This is not true for 
air showers parallel to the geomagnetic field. In such cases, systematic deviations of 
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Figure 7: Polarization vectors of the transverse current and the net charge excess contributions in 
the plane perpendicular to the shower axis. Due to the fact that the polarization vectors of these two 
contributions are not always oriented in the same direction, their combination can be constructive or 
destructive in function of the antenna position. For the east-west polarization of the electric field (as it is 
given in Fig|4] and Figjs] left) , the signal amphtudes observed by antennas located on the east side of the 
plane containing the ground shower core and the geomagnetic field appear finally higher than the signal 
amplitudes of antennas located on the other side. 
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Figure 8: Left: Antennas positions used for ground footprint calculation. Center and right: Ground 
footprint of the energy deposited by a vertical 10"'^^ eV air shower below 85 MHz, for the two horizontal 
polarizations. The east-west polarization shows a slight east-west asymmetry which is due to the net 
charge excess contribution (see text for more details). 
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particles due to the geomagnetic field disappear and the transverse current contribution 
vanishes. To understand the behavior of the two dominant contributions combination, 
we show in FigJT] the polarization vectors of the transverse current and the net charge 
excess contributions in the plane perpendicular to the shower axis. Due to the fact that 
the polarization vectors of these two contributions are not always oriented in the same 
direction, their combination can be constructive or destructive. In the configuration of the 
Auger site for instance, the summation is then constructive for the case of antennas located 
on the east side of the plane containing the ground shower core and the geomagnetic field. 
On the west side of this plane the summation is destructive. This asymmetry effect is 
also observable in an other form in FigjS] where we present the ground footprint of the 
two horizontal polarizations. FigjS] are obtained using a dense array composed of 145 
antennas (cf Figls] left) for which we calculate the energy deposited below 85 MHz, for 
the two horizontal polarizations. The results obtained here confirm results previously 
obtained and discussed in [TTf [H]. 
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Figure 9: From [517. Top: Comparison of the western polarization unlimited bandwidth radio pulses 
for vertical showers with a primary energy of 10^^ eV for REAS3 (left) and MGMR (right). Bottom: 
Comparison of the absolute field strength unlimited bandwidth pulse amplitude lateral distribution for 
vertical showers with a primary energy of 10"'^^ eV for REAS3 (left) and MGMR (right). These figures 
can be directly compared with Fig|4]and Fig|5] 



4.4. Comparison between SELFAS2, REAS3 and MGMR 

In a recent paper [20j, a direct comparison is done between the two main approaches 
discussed in the section 2. They show in this paper that after great disagreements, the 
two models finally give results which converge in the same direction. We extracted some 
figures of this paper to make a direct comparison with the results obtained with SELFAS2. 
In Figj9] we show the results obtained by REAS3 and MGMR for the same air shower 
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simulated with SELFAS2 in subsection 4.1. These figures can be directly compared with 
Figj4] and FigjS} The agreement between the three models is satisfactory but it still 
remains some differences which can be due to the electric field formalism adopted and 
to the shower description. For a better comparison, the simulations performed with 
REAS3 and MGMR are based on EAS characteristics coming from the same CORSIKA 
shower output file. SELFAS2 doesn't use any simulated shower, but only the analytical 
geometrical distributions previously described (based on CORSIKA simulation). The 
agreement between the three models shows that their implementation in SELFAS2 gives 
promising results. 

5. Conclusion 

With SELFAS2, we successfully managed to construct a complete autonomous tool 
able to help detailed understanding in the mechanisms involved in the air showers radio 
emission. The concept of universality of air showers feature permits to avoid the heavy 
use of air shower generators and to obtain acceptable results after few minutes for one 
antenna simulated. Furthermore, we confirmed that EAS radio emission is mainly due to 
two different contributions: the time variation of the transverse current (corresponding 
to the third term in EqjT7|) and the time variation of the EAS macroscopic charge excess 



(corresponding to the second term in Eq, 17). The polarization patterns of these contri 



butions are different leading to asymmetries in both the EAS ground footprint and the 
lateral profile. It would be very interesting to search for such asymmetries in the data 
(CODALEMA, LOPES, RAuger) in order to make progress in the understanding of the 
radio emission. The quantification of the experimental radio asymmetry should be a good 
estimator of the electrons-positrons asymmetry. 

Today, we can clearly say that air shower radio simulations converge in the same 
direction, showing good agreements even using different approaches. Coupled analysis 
between data and SELFAS2 are currently in progress to test the validity of theoretical 
predictions. 
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